(thermodynamic  properties 

OF  CESIUM  UP  TO  1500°  k1 


by  Sheldon  Heimel 

Lewis  Research  Center 
Cleveland,  Ohio 


reproduced  from 
best  avaiiable  copy 

20010906  073 


NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 


WASHINGTON,  D.  C. 


JULY  1965 


NASA  TN  D-2906 


THERMODYNAMIC  PROPERTIES  OF  CESIUM  UP  TO  1500°  K 

By  Sheldon  Heim  el 

Lewis  Research  Center 
Cleveland,  Ohio 


NATIONAL  AERONAUTICS  AND  SPACE  ADMINISTRATION 


For  sale  by  the  Clearinghouse  for  Federal  Scientific  and  Technical  Information 
Springfield,  Virginia  22151  -  Price  $2.00 


CONTENTS 


Page 

SUMMARY .  1 

INTRODUCTION .  2 

PRESENTATION  OF  TABLES .  3 

Assigned  Reference  State .  ^ 

Assigned  Enthalpy  Values  .  ^ 

Heats  of  Formation .  3 

Equilibrium  Constants . 

SELECTION  OF  INITIAL  DATA .  10 

Crystal .  10 

Liquid .  H 

Monomer  (Gas) .  12 

Dimer  (Gas) .  13 

Vapor  Pressures .  13 

Heat  of  Dissociation  at  0®  . .  I*! 

Heat  of  Condensation  of  Monomer  at  0®  K .  14 

THERMODYNAMIC  ANALYSIS .  14 

Derivation  of  (AH^)^ .  1^ 

Partial  Pressure  of  Monomer .  13 

RESULTS  OF  CALCULATING  (AH^)^ .  10 

CALCULATION  OF  EQUILIBRIUM  MIXTURE  PROPERTIES  ON  SATURATION 

LINE .  10 

Pressures  and  Weight  Fractions  of  Gaseous  Species .  18 

Enthalpy  of  Vaporization  and  Entropy  of  Vaporization  .  18 

Enthalpy,  Specific  Heat,  and  Entropy  of  Cesium  Vapor .  18 

Molecular  Weight  of  Mixture .  20 

Specific  Volume  and  Density  of  Mixture .  21 

DISCUSSION  OF  VAPOR  PRESSURES  . 21 

Vapor  Pressures  Above  Liquid  .  21 

Vapor  Pressures  Above  Crystal .  22 

APPENDIX  -  SYMBOLS  .  25 

REFERENCES .  27 


THERMODYNAMIC  PROPERTIES  OF  CESIUM  UP  TO  1500°  K 

by  Sheldon  Helmel 
Lewis  Research  Center 

SUMMARY 


persistent  tables  of  thermodynamic  propertlejg  of  elemental  cesium  were  compiled 
for  both  the  pure  species  and  the  equilibrium  vapor  on  the  saturation  line  using  selected 
values  of  10  500  calories  per  mole  for  the  heat  of  dissociation  of  dimer  and  -18  920  cal¬ 
ories  per  mole  for  the  heat  of  condensation  of  monomer.  The  equilibrium  vapor  prop¬ 
erties  and  properties  of  condensed  cesium  (Cs)  are  given  to  1500°  K,  while  properties  of 
gaseous  monomer  and  dimer  are  given  to  2500°  K.  f,  ^ 

The  table  of  mixture  properties  on  the  saturation  line  gives  the  following  thermo¬ 
dynamic  functions:  temperature  T,  total  vapor  pressure  P,  weight  fraction  of  dimer  in 
vapor  X2,  molecular  weight  Mg,  specific  volume  v,  density  p,  enthalpy  of  condensed 
phase  he,  enthalpy  of  vaporization  Ah^^p,  enthalpy  of  vapor  hg,  entropy  of  condensed 
phase  Sg,  entropy  of  vaporization  As^g^p,  entropy  of  vapor  Sg,  and  heat  capacities 
(equilibrium  and  frozen)  (cp)gq  and  (cp)j^.  Assigned  enthalpies  are  relative  to  enthalpy 

of  the  crystal  at  298. 15°  K,  which  is  taken  to  be  zero.  The  tables  of  properties  for  the 
pure  species  give  the  following  thermodynamic  functions  for  the  standard  state:  heat 
capacity  at  constant  pressure  C°,  sensible  enthalpy  -  H^,  entropy  S^,  sensible 
free  energy  -  Hq,  the  sum  of  sensible  enthalpy  at  T°  K  and  chemical  energy  at 
0°  K  H^,  and  values  of  enthalpy  changes  and  logarithms  of  the  equilibrium  constants. 
The  latter  two  functions  are  given  for  the  reactions  of  formation  of  Cs(c),  CSj^(g),  and 
Cs2(g)  from  the  assigned  reference  state  Cs(c)  ((AH^)^  and  logjQ  Kj)  and  from  Csj(g) 

(AH^  and  log^Q  K). 

/Thermodynamic  functions  for  the  gases  were  generated  from  atomic  and  molecular 
data,  whereas  the  functions  of  the  condensed  phase  were  based  on  selected  experimental 
data  that  were  smoothed  and  made  self-consistent.\  . 


TABLE  I.  -  THERMODYNAMIC  PROPERTIES 


OF  CESIUM  (CRYSTAL)  UP  TO  100°  K 
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95 
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tables  III  (Csj)  and  IV  (CS2),  data  are 
tabulated  at  100®  intervals  from  0®  to 
2500®  K  for  the  ideal  gases.  Table  V 
(p.  8)  contains  the  molecular  con¬ 
stants  used  to  calculate  thermody¬ 
namic  functions  for  dimeric  cesium. 

Table  VI  (p.  8)  contains  recom¬ 
mended  boiling  points  of  liquid 
cesium  calculated  from  the  final  set 
of  thermodynamic  properties.  These 
boiling  points  correspond  to  pres- 

_  O 

sures  from  10”  to  10  atmospheres 
at  every  power  of  10.  For  compar¬ 
ison,  boiling  points  obtained  from 
equation  (33)  are  also  given. 

Table  VII  (p.  9)  gives  the  follow¬ 
ing  thermodynamic  properties  of  the 
equilibrium  gaseous  mixture  on  the 
saturation  line:  T,  P,  pj,  P2, 

Mg,  V,  p,  h^,  ^h^ap’  ®c’  ^®vap’ 

s  ,  (c  )  ,  and  (c  )  .  The  values  of 

^  P  eq  P  fr 

P2>  P  given  in 

floating-point  notation,  where  the 
decimal  number  is  to  be  multiplied  by 


10  raised  to  the  power  of  the  sign  and  the  two  digits  following  the  letter  E;  for  example, 
0. lOOOOOE-03  is  0. 100000x10”^  =  0. 0001. 


Assigne(j  Reference  State 

The  assigned  reference  state  is  Cs  (crystal,  liquid).  The  crystal  is  the  reference 
state  up  to  the  melting  point,  and  the  liquid  phase  is  the  reference  state  above  the  melt¬ 
ing  point. 


Assigned  Enthalpy  Values  Hj 

For  some  applications  (ref.  1),  it  is  convenient  to  combine  sensible  enthalpy 
-  Hq  and  chemical  energy  Hq  into  one  numerical  value,  H^.  The  arbitrary  base 
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The  reference  element  is  crystal  cesium  up  to  301. 8^  K  and  liquid  cesium  above  301. 8^ 
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TABLE  V.  -  MOLECULAR 


TABLE  VI.  -  BOILING  POINTS  OF  LIQUID 


CONSTANTS  FOR  DI¬ 
ATOMIC  CESIUM 


CESIUM  FROM  10“°  TO  10  ATMOSPHERES 


Molecular  weight 
Symmetry  number 
Electronic  state^ 

Statistical  weight 
-1^ 

WpYp,  cm'-^ 

®  ®  lb 

Bp,  cm-l 
®  -lb 

a„,  cm  ^ 

®  lb 

Dg,  cm-1 


Constant  from  ref. 
’constant  from  ref. 


265.82 

2 

1 

41.990 
0.  08005 
-0. 0001643 
0.01272 
0.000035 
0.467X10'® 


^atm 

Boiling  points, 

Recommended 

From  equation  (33)^ 

-8 

317.  8 

317.8 

-7 

345.  8 

345.8 

-6 

379.  2 

379.4 

-5 

420.  1 

420.3 

-4 

471.  2 

471.4 

-3 

537.0 

537.  2 

-2 

625.0 

625.0 

-1 

749. 1 

748.  7 

0 

937.  7 

937.  3 

1 

1264.0 

1264. 5 

■10  ^atm 


^-4053.30^7  Q4453  _  q  915282  logjQ  T 


for  assigning  values  to  the  enthalpy  of  Cs(c)  was  a  value  of  zero  at  298. 15°  K.  Since 

table  n  gives  (H298  15  -  Hq)  =  1807.  2  calories  per  mole,  then  (H^) 

Cs(c)  Cs(c) 

-1807.  2  calories  per  mole.  From  table  II,  (AH^)  =  (H^)  -  (H^)  = 

c  ^  Cs(c)  ^  Csj(g) 

-18  920.  0  calories  per  mole.  Therefore,  (H^)  =  -1807.  2  -  (-18  920.  0)  = 

^  Csj(g) 

17  112.  8  calories  per  mole.  From  table  IV,  (AH^)  =  -D^  =  (H^)  -  2(H^) 

°2  O  ®Cs2(g)  ^Csi(g) 

-10  500  calories  per  mole.  Therefore, 

(H^)  =  -10  500.  0  +  2(17  112.  8)  =  23  725.  6  cal/mole 

^  Cs2(g) 


Heats  of  Formation 

Two  sets  of  values  for  heats  of  formation  are  given  in  tables  II  to  IV.  The  first  set 
(col.  7)  is  for  the  formation  of  the  given  species  from  Cs(c),  (AH^)^;  the  second  set 

(col.  9)  is  for  the  formation  from  Csj(g),  (AH^).  For  Cs2(g)  at  298. 15°  K,  for  ex¬ 
ample. 


298. 15 


)  =(H° 


298. 15' 


Cs2(g) 


^298. 15' 


Cs(c) 


8 


TABLE  VH.  -  THERMODYNAMIC  PROPERTIES  OF  CESIUM  ON  SATURATION  LINE 
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16>2  -  <H^98.  15>c,^(g) 


(2) 


Equilibrium  Constants 


Two  sets  of  logarithms  of  the  equilibrium  constants  for  the  two  formation  reactions 
discussed  previously  are  also  listed  in  tables  II  to  IV.  The  equilibrium  constant  Kj  for 
formation  from  the  assigned  reference  element  (col.  8)  is  obtained  from  the  standard 
free  energy  change  (AF^)  by  means  of  the  equation 


(AF°) 

log.oKf= - - -  (3) 

^  2.  3025851  RT 

The  equilibrium  constant  K  for  formation  from  the  atomic  gas  is  obtained  from  a  sim¬ 
ilar  equation 


AF^ 

log  K  - - ^ -  (4) 

2.  3025851  RT 

Atomic  weights,  the  universal  gas  constant  and  the  constants  used  in  the  evaluation 
of  the  entropy  constant  were  the  same  as  those  used  in  reference  1. 

SELECTION  OF  INITIAL  DATA 
Crystal 

Heat  capacities  of  the  crystal  (table  I,  p.  4)  were  derived  by  smoothing  the  experi¬ 
mental  data  from  the  sources  given  in  table  VIII. 

The  following  procedure  was  used  to  obtain  smooth  data  from  0°  to  100°  K.  A 
least-squares  fit  of  the  data  of  reference  2  was  used  in  the  temperature  range  from  0°  to 
4°  K.  Data  from  4°  to  20°  K  were  obtained  from  a  curve  drawn  through  the  data  of  ref¬ 
erence  2  at  4°  K  and  through  the  data  of  references  3  and  4  at  20°  K.  Data  from  20°  to 
100°  K  were  obtained  from  a  least-squares  fit  of  the  data  of  reference  4  in  that  tempera¬ 
ture  range. 

Between  100°  and  210°  K,  the  data  of  reference  4  are  approximately  linear  in  tem- 
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TABLE  Vm.  -  REFERENCES  FOR  HEAT 
CAPACITIES  OF  CRYSTAL 


Temperature 

range, 

°K 

Source 

0.  3  to  4 

Ref.  2 

8  to  12 

Private  communication  from 

Professor  D.  C.  McCollum 

of  University  of  California 

in  Riverside,  California 

and  ref.  3 

20  to  210 

Ref.  4 

perature.  The  heat  capacities  at  lOO'^  and 
210°  K  yield  the  equation 


Cp  =  5.  76  +  0.  004  T 


(5) 


Above  210°  K,  thermal  effects  associated  with 
the  cesium-oxygen  system  cause  anomalous  in¬ 
creases  in  heat  capacity  and  therefore  these 
data  were  not  used.  Equation  (5)  was  used  to 
extrapolate  heat  capacities  up  to  the  melting 
point.  The  melting  point  of  301.  8°  K  and  the 
heat  of  fusion  of  520. 1  calories  per  mole  were  taken  from  reference  5. 

Reference  4  gives  a  melting  point  of  300.  5°  K  and  the  total  heat  absorption  between 
the  crystal  at  262.  5°  K  and  the  liquid  at  300.  5°  K.  If  this  heat  is  compared  with  the 
sensible  enthalpy  change  of  the  crystal  between  these  two  temperatures,  as  calculated 
from  equation  (5),  the  implied  heat  of  melting  is  514  calories  per  mole.  This  value  is 
in  reasonable  agreement  with  the  value  adopted. 

The  enthalpy  from  the  linear  equation  for  C°,  when  combined  with  the  heat  of  fusion, 
gave 


^%01. 8^cs(liquid) 


(H° 


298. 15 


Cs(crystal) 


=  545.  5 


cal/  mole 


(6) 


Liquid 

A  number  of  empirical  equations  have  been  used  to  satisfactorily  fit  enthalpy  data  as 
a  function  of  temperature  for  various  substances;  however,  heat  capacity  data  derived 
from  differentiation  of  the  enthalpy  equations  are  often  very  xmreliable.  For  any  par¬ 
ticular  species,  various  empirical  enthalpy  equations  that  might  be  chosen  yield  con¬ 
siderably  different  sets  of  heat  capacities. 

Two  sets  of  experimental  enthalpy  data  were  considered:  that  of  references  6  and  7. 
The  curve  of  reference  6  exhibits  an  anomalous  bump  from  351°  to  620°  K,  followed  by 
a  straight  line  to  1238°  K.  It  is  recommended  in  reference  6  that  just  the  straight-line 
portion  be  used,  which  has  a  standard  deviation  a  =  1.  5  percent.  For  the  cubic  equa¬ 
tion  (from  340°  to  1176°  K)  of  reference  7,  a  =  1. 1  percent. 

The  implied  C°  =  7.  25  calories  per  mole  per  °K  of  reference  6  differs  greatly 

£r 

from  the  heat  capacities  derived  from  the  enthalpy  curve  of  reference  7  and  listed  in 
table  K. 

According  to  reference  7,  the  curves  of  heat  capacity  against  temperature  of  liquid 
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TABLE  IX.  -  HEAT  CAPACITIES  cesium  and  other  alkali  metals  have  a  parabolic  shape; 


DERIVED  FROM  ENTHALPY 
CURVE  OF  REFERENCE  7 


Absolute 

temperature, 

T, 

Heat  capacity  at 
constant  pressure  for 
standard  state, 

cal/(moIe)(°K) 

301.8 

10.35 

500 

7.  28 

700 

6.23 

1000 

8.53 

1200 

12.64 

1500 

22.67 

that  is,  heat  capacity  initially  decreases  with  increas¬ 
ing  temperature,  reaches  a  minimum,  and  then  in¬ 
creases  continuously.  However,  a  parabolic  shape  for 
heat  capacity  of  these  species  is  only  the  result  of 
selecting  a  cubic  to  represent  experimental  enthalpy 
data. 

It  was  thought  advisable  to  use  a  straight  line  for 
enthalpy  because  (1)  the  straight  line,  in  the  standard 
deviation  sense,  represents  cesium  enthalpies  about 
as  well  as  the  cubic  does  and  (2)  the  cubic  has  the  dis¬ 
advantage  of  implying  very  large  heat  capacities  for 
even  a  few  hundred  degrees  extrapolation. 

The  actual  straight  line  used  in  this  report  was 
based  on  the  straight-line  data  of  reference  6  but  was 


constrained  to  go  through  the  previously  calculated  relative  heat  content  of  545.  5  cal¬ 


ories  per  mole  at  the  melting  point  (eq.  (6)).  This  gives  the  following  equation: 


(H°)  , 

^  Cs(liquid) 


^^298  15^ 

298.  10  cs{crystal) 


=  -1704.  5  +  7.455  T 


(7) 


Equation  (7)  implies  a  constant  value  of  C°  =  7.  455  calories  per  mole  per  °K, 

Had  the  data  of  reference  7  been  fitted  to  a  straight  line  and  similarly  constrained  to 
go  through  the  heat  content  of  545.  5  calories  per  mole,  a  constant  value  of  C°  =  7.  595 
calories  per  mole  per  °K  would  have  resulted.  The  enthalpies  would  be,  on  the  average, 
only  about  1  percent  higher  than  those  derived  from  the  data  of  reference  6. 

The  enthalpy  data  of  reference  8  are  about  6  to  7  percent  lower  than  the  data  of  both 
references  6  and  7  over  the  range  from  about  600*^  to  1300^  K;  therefore,  they  were  not 
used. 


Monomer  (Gas) 

The  partition  function  used  to  compute  the  thermodynamic  properties  of  Csj  was 
truncated  by  the  temperature -dependent  cutoff  technique  used  in  reference  1.  The  re¬ 
sults  of  these  computations  are  exactly  the  same  up  to  2200^  K  as  would  be  obtained  by 
simply  using  all  levels  given  in  reference  9  with  no  cutoff.  Above  2200*^  K  enthalpy  and 
entropy  differ  only  in  the  last  place  computed  and  heat  capacity  in  the  last  or  second  last 
place. 
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Dimer  (Gas) 


The  method  used  to  compute  the  thermodynamic  properties  of  Cs2  is  the  method 
for  diatomic  molecules  described  in  reference  1.  The  constants  for  the  dimer  are 
shown  in  table  V  (p.  8).  They  were  taken  from  references  10  and  11. 


Vapor  Pressures 

During  the  period  1913  to  1937,  the  investigations  reported  in  references  12  to  18 
obtained  vapor  pressure  data  for  cesium  in  the  temperature  range  from  238°  to  670°  K. 
This  early  work  was  evaluated  in  reference  19  where  it  was  concluded  that,  except  for 
the  results  of  references  17  and  18,  the  early  work  is  not  too  reliable. 

The  vapor  pressures  in  the  455°  to  589°  K  range  (ref.  17)  were  used  in  reference  19 
to  generate  a  vapor  pressure  equation.  Inasmuch  as  the  data  of  reference  17  are  only 
relative  pressure  data  (obtained  by  the  magneto -optical  method)  and  are  converted  to 
absolute  pressure  data  by  relating  them  to  the  unreliable  data  of  reference  13,  it  seems 
prudent  to  disregard  the  data  of  reference  17. 

Recently,  the  first  experimental  measurements  in  the  temperature  range  that  in¬ 
cluded  the  boiling  point  (742°  to  1199°  K)  were  reported  in  reference  20.  These  data 
check  reasonably  well  with  the  later  data  of  reference  6  (727°  to  1334°  K). 

In  view  of  the  previous  observations,  the  data  of  references  18  (vapor  pressures  in 
the  liquid  range),  20,  and  6  were  combined  and  fitted  by  the  least-squares  technique. 

All  the  data  were  given  the  same  weight.  The  results  of  several  empirical  equations 
were  compared.  The  following  equation,  which  is  in  the  form  of  the  Kirchhoff  equation, 
was  accepted  because  it  gave  as  satisfactory  results  as  the  other  forms  considered: 

logjo  P  =  +  5.  71342  -  0.  519781  log^Q  T  (8) 

If  only  vapor  pressures  were  considered  in  this  report,  then  the  smoothed  experi¬ 
mental  vapor  pressures  obtained  from  equation  (8)  would  be  accepted.  As  pointed  out  in 
the  INTRODUCTION,  however,  these  smoothed  data,  when  combined  with  other  thermo¬ 
dynamic  data,  lead  to  inconsistencies  in  the  heat  of  condensation  of  the  monomer  at  0°  K 
(AH^)  .  Therefore,  the  data  from  equation  (8)  are  used  only  for  preliminary  analysis, 

whereas  the  final  recommended  vapor  pressures  are  calculated  as  explained  in  detail  in 
the  section  Pressures  and  Weight  Fractions  of  Gaseous  Species. 
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Heat  of  Dissociation  at  0°  K 


An  upper  bound  (4020  cm”^)  for  the  heat  of  dissociation  at  0°  K(Dq)  is  estimated  in 
reference  21  by  analytical  extrapolation  of  the  ground  state  vibrational  intervals.  A 
lower  bound  (3318  cm'^)is  estimated  in  reference  21  by  making  certain  assumptions  as  to 
the  final  states  of  the  excited  atoms  produced.  The  average  of  these  two  values,  which  is 
equivalent  to  10  500±1000  calories  per  mole  of  dimer,  is  recommended  in  this  reference. 

The  estimate  of  0.  453  eV  or  10  450  calories  per  mole  from  the  correlation  of  force 
constants,  equilibrium  internuclear  distances,  and  heats  of  dissociation  of  the  alkali 
metals  is  obtained  in  reference  22.  A  technique  for  deducing  a  value  of  Dq  where  Dq 
was  not  well  known  is  suggested  in  reference  23.  This  technique  gave  an  optimum 
Dq  =  10  000  calories  per  mole  which  is  in  fairly  close  agreement  with  the  spectroscopic 
value  of  10  500  calories  per  mole. 

The  calculation  of  (AH^)  ,  which  is  discussed  in  the  next  section,  depends  in  part 

on  the  value  of  Dq,  While  no  value  of  Dq  that  was  tried  gave  a  constant  value  of 
(AHq)  over  the  entire  temperature  range  of  301.  8*^  to  1500®  K,  values  of  Dq  from 

about  10  000  to  11  500  calories  per  mole  gave  the  most  nearly  constant  value  of  (AH^)  . 

Based  on  this  analysis,  the  conclusion  was  that  the  spectroscopic  value  was  ap-  ® 
proximately  correct,  and  the  value  Dq  =  10  500  calories  per  mole  was  selected. 

Heat  of  Condensation  of  Monomer  at  0°  K 

The  values  of  (Ab£)  in  calories  per  mole  are  given  in  a  number  of  publications  as 
c 

follows:  -19  048  (ref.  24),  -18  790  (ref.  19),  -19  050  (ref.  25),  and -19  035  (ref.  26). 
The  aforementioned  values  of  (AH^)  are  based  on  the  pre-1962  vapor  pressures. 

Therefore,  this  report  will  recommend  a  "best”  (AH^)  of  -18  920  on  the  basis  of  the 
newer  vapor  pressures,  as  discussed  in  the  following  section. 


THERMODYNAMIC  ANALYSIS 


As  discussed  in  the  INTRODUCTION,  the  procedure  for  the  selection  of  a  value  for 

(AH^  is  to  find  that  (AH^)  for  which  the  difference  between  the  calculated  and 
c  c 

smoothed  experimental  vapor  pressures  is  a  minimum.  As  a  first  step,  it  is  necessary 
to  compute  (AH^)  for  each  temperature  in  the  range  of  interest. 
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Derivation  of  (AHn) 

c 

The  zero-degree  heat  of  condensation  of  monomer  was  computed  by  means  of  the 
relation 

(AH°)  =  (AF°)  -  A(f2,  -  H°)  (9) 

c  c  c 

where 

(AF°)  =(F°)  -(F°)  =RTlnpj  (10) 

^  c  ^  Cs(c)  ^  Csj(g) 

and 

A(F°  -  Hg)  =  (F°  -  Hg)  -  (F°  -  H^)  (11) 

Partial  Pressure  of  Monomer 

The  partial  pressure  of  monomer  pj  in  atmospheres  is  obtained  from  the  simul¬ 
taneous  solution  of  Dalton’s  Law  and  the  equilibrium- constant  equation  for  the  dimeriza¬ 
tion  reaction 

2  CSj  =  CSg 

The  equations  are 

P  =  Pj  +  P2  (12) 

and 

K2  =  ^  (13) 

^  2 
Pi 

where  P  is  the  total  pressure,  P2  is  the  partial  pressure  of  dimer,  and  K2  is  the 
equilibrium  constant  of  dimerization.  Values  of  K2  are  obtained  from  table  IV  (p.  7). 
The  solution  of  equations  (12)  and  (13)  is 

15 


-1 + yiPK^TT 


(14) 


Pi  = 


2K2 


Therefore,  apparent  values  of  (AH^)  can  be  computed  at  every  temperature  of  interest 

■  c 

by  starting  with  given  values  of  Dq,  values  of  for  monomer,  dimer,  and 

condensed  phase,  and  smoothed  experimental  vapor  pressures  Ppjj*  This  procedure 
can  be  summarized  as  follows: 


(fJ  -  H°) 


Cs(c) 
(Table  II) 


(F°  -  H^) 


Cs^(g) 
(Table  IH) 


Equation  (11) 


A(F®  -  hJ) 


^2 

(Table  IV) 


Equation  (9) 


(A<) 


^FIT 
(Eq.  (8)) 


Equation  ( 14) 


Pi 


RESULTS  OF  CALCULATING  (AHR) 

^  c 

The  apparent  (AHq)  values  computed  by  means  of  this  procedure  from  equation  (9) 
^  c 

are  given  in  table  X.  The  maximum  difference  is  about  1  percent. 

If  some  single  arbitrary  value  of  (AH^)  is  selected,  then  for  each  of  the  previous 

temperatures  there  will  be  a  difference  between  the  smoothed  experimental  vapor  pres¬ 
sure  PjTfp  and  the  vapor  pressure  Pq^LC  with  this  arbitrary  (AHq) 

value.  ^ 

For  this  analysis,  an  error  function  may  be  defined  as  [Ppyj-  -  PcALcl^^FIT 
averaged  over  the  previous  number  of  temperatures.  The  procedure  for  selecting  a  value 
of  (AHo)(.  is  to  find  that  value  of  (AHq)^  which  gives  the  minimum  error  function. 

The  results  of  calculating  the  error  function  for  various  values  of  (AHo)^  are  given 
in  table  XI.  The  optimum  (AHq)^  was  taken  to  be  -18  920  calories  per  mole  to  the 
nearest  5  calories. 

The  fact  that  there  is  still  about  4  percent  error  between  the  smoothed  experimental 

vapor  pressure  and  vapor  pressure  calculated  with  the  optimum  (AHq)  ,  indicates  that 

c 
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TABLE  X.  -  APPARENT  VALUES  OF  -(AH^) 


c 


AS  A  FUNCTION  OF  TEMPERATURE 


Tempera¬ 

ture, 

T, 

°K 

cal/mole 

Tempera¬ 

ture, 

T, 

°K 

-(AH^)  , 
^  c 

cal/mole 

301.8 

18  791.  6 

1000 

18  965. 1 

400 

18  872.5 

1100 

18  954.  3 

500 

18  921.7 

1200 

18  937.0 

600 

18  949.  3 

1300 

18  913. 1 

700 

18  963.  6 

1400 

18  881.9 

800 

18  969.  7 

1500 

18  843.4 

900 

18  970.  0 

TABLE  XI.  -  ERROR  FUNCTION 
FOR  VARIOUS  VALUES 


OF  -(AH®) 

c 


-{AH°K 
cal/  mole 

Error  function, 
|ap  |/p 

18  850 

0, 0538 

18  900 

.0397 

18  910 

.0381 

18  920 

.0370 

18  925 

.0370 

18  930 

.0372 

18  950 

.0389 

18  970 

.0458 

discrepancies  still  remain  which  may  be  in  any  of  the  following  data: 

(a)  Experimental  vapor  pressures 

WdS 

(c)  Molecular  constants  of  dimer  (and,  consequently,  its  thermodynamic  properties) 

(d)  Heat  content  (and  other  derived  thermodynamic  properties)  of  the  condensed  phase 

(e)  (AH°)^ 

The  4  percent  error  may  also  be  due,  in  part,  to  the  presence  of  a  higher  polymer,  for 
example,  trimer  or  tetramer;  however,  this  possibility  has  not  been  considered  in  this 
report. 


CALCUUMION  OF  EQUILIBRIUM  MIXTURE  PROPERTIES 
ON  SATURATION  LINE 

As  explained  in  the  INTRODUCTION,  selection  of  an  optimum  value  of  (AHq)  makes 

c 

it  possible  to  compute  a  set  of  vapor  pressures  consistent  with  the  sensible  free  energies 
of  the  pertinent  species.  These  vapor  pressures  and  other  equilibrium  properties  on  the 
saturation  line  (e.  g. ,  enthalpy,  specific  heat,  and  entropy)  are  given  in  table  Vn  (p.  9). 
The  equations  used  to  obtain  these  properties  are  given  in  the  following  sections. 


Pressures  and  Weight  Fractions  of  Gaseous  Species 


By  starting  with  the  optimized  (AH^)  and  the  sensible  free  energies  of  monomer 

c 

and  condensed  phase,  monomer  pressure  was  computed  with  equations  (9)  to  (11). 
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Dimer  pressure  was  computed  from  equation  (13).  The  total  pressure 
obtained  from  the  sum  of  Pj^  and  pg^ 

The  weight  fraction  of  monomer  was  obtained  from  the  mole  fraction  of 
monomer  Nj^;  thus, 

Pi 

N,  =-^  (15) 

^  P 


N. 


Xi  = 


2  -  N. 


(16) 


The  weight  fraction  of  the  dimer  X2  was  obtained  from  the  defining  equation 

Xi  +  X2  =  1 


(17) 


Enthalpy  of  Vaporization  and  Entropy  of  Vaporization 

The  heat  of  vaporization  in  calories  per  gram  from  the  condensed  state  to  the 
equilibrium  mixture  on  the  saturation  line  was  obtained  from 


ap 


2x^(AHt)j^  ^  +  X2(AHT)f^  2 


(18) 


M, 


where  M2  is  the  molecular  weight  of  dimer  (265.  82  g/mole).  The  heats  of  formation 
are  obtained  from  tables  III  and  IV  (pp.  6  and  7). 

The  entropy  of  vaporization  in  calories  per  gram  is  obtained  from  equation  (18)  by 
means  of 


As 


vap 


T 


(19) 


Enthalpy,  Specific  Heat,  and  Entropy  of  Cesium  Vapor 


The  enthalpy 


‘g 


in  calories  per  gram  of  an  equilibrium  mixture  can  be  defined  by 
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(20) 


Equation  (20)  for  saturation  conditions  gives  the  same  results  as  the  usual  equation  for 
the  enthalpy  of  a  mixture 


where  h  is  the  enthalpy  of  the  condensed  phase  in  calories  per  gram  (table  VII,  p.  9). 

^  o 

The  equilibrium  specific  heat  at  constant  pressure  in  calories  per  gram  per  K  is 

defined  by 


Using  equation  (20)  results  in  equation  (21)  becoming 


or 


From  equation  (17),  it  follows  that 


Therefore, 


(c 


(hg  -  hj)  = 


'9x, 


8T 


(AH^) 

^  2 

Mg 


(21) 


(22) 


(23) 


(24) 


(25) 
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Values  of  (AH^)  are  tabulated  in  table  IV  (p.  7). 

T  2 

calculated  from  the  following  equation: 


Values  of  (0X2/9T)p  may  be 


^2(2  -  X2){1  -  X2) 

Rt2  2 


(26) 


Equation  (26)  may  be  derived  from  the  following  expression,  which  is  equivalent  to 
equation  (13): 


X2(2  -  X2) 

4P(1  -  x/ 


and  from  the  van't  Hoff  isobar,  namely, 


(27) 


(28) 


The  entropy  of  the  gas  mixture  may  be  obtained  from 


Sg  =  ^Svap  + 


(29) 


where  s^  is  the  entropy  of  the  condensed  phase  in  calories  per  gram  per  °K  (table  Vn, 
p.  9)  and  As^^p  is  obtained  from  equation  (19).  Equation  (29)  gives  the  identical 
results  as  the  usual  equation  for  the  entropy  of  a  mixture. 


Molecular  Weight  of  Mixture 

An  expression  for  the  molecular  weight  of  the  mixture  Mg  (in  g/ mole)  in  terms  of 
the  molecular  weight  and  weight  fraction  of  the  dimer  may  be  obtained  as  follows: 


1 

ni  +  n2 


M2 

2  -  X2 


(30) 
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where  n,  and  Ug  are  the  number  of  moles  of  mixture,  monomer,  and  dimer  per 
gram  of  mixture,  respectively. 


Specific  Volume  and  Density  of  Mixture 

The  specific  volume  v  in  cubic  centimeters  per  gram  may  be  obtained  from  the 
ideal  gas  law 


_  RT  _  82.  05971  T 

PM„  PM„ 

g  g 

The  density  p  in  grams  per  1000  cubic  centimeters  is,  therefore, 

V 


(31) 


(32) 


DISCUSSION  OF  VAPOR  PRESSURES 
Vapor  Pressures  Above  Liquid 

-8 

Table  VI  contains  boiling  points  of  liquid  cesium  at  saturation  pressures  from  10" 
to  10  atmospheres  at  every  power  of  10.  These  data  were  calculated  in  two  ways: 

(1)  From  the  recommended  thermodynamic  data,  as  explained  in  the  section 

THERMODYNAMIC  ANALYSIS 

(2)  By  a  least-squares  equation  fitted  to  the  recommended  P^^l^  values  with  a 

standard  deviation  of  0.  5  percent 

Patm  =  +  ’7-  04453  -  0.  915282  logjQ  T  (33) 

On  the  average  the  two  sets  of  boiling  points  differ  by  about  0.  2°  K.  It  should  be 
emphasized  that  the  vapor  pressure  equation  (eq.  (33))  is  recommended  rather  than 
equation  (8). 

Figures  1  and  2  compare  the  experimental  vapor  pressures  of  liquid  cesium  with  a 
curve  drawn  through  the  recommended  vapor  pressures.  The  agreement  of  the  two  sets 
is  seen  to,  be  very  good.  Figure  1  covers  the  temperature  range  from  303*^  to  1334®  K. 
Figure  2  covers  only  the  range  from  729®  to  1334®  K  to  facilitate  detailed  comparison. 
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Patm 


Figure  1.  -  Comparison  of  experimental  and  calculated  vapor  pressures  above  liquid  cesium  (303°  to  1334°  K). 


Figure  3.  -  Comparison  of  experimental  and  calculated  vapor  pressures  above  crystalline  cesium. 

Subsequent  to  the  completion  of  the  analysis,  a  vapor  pressure  equation  in  the  tem¬ 
perature  range  928°  to  1558°  K  was  presented  in  reference  27.  On  the  average,  these 
vapor  pressures  are  about  2  percent  lower  than  those  recommended  in  this  report. 


Vapor  Pressures  Above  Crystal 

Experimental  vapor  pressures  above  the  crystal  are  available  from  reference  18  at 
temperatures  from  238°  to  298°  K.  In  order  to  be  consistent  with  the  thermodynamic 
data  adopted  for  this  report,  however,  vapor  pressures  from  100°  to  301.  8°  K  were  cal¬ 
culated  in  the  same  manner  as  were  the  vapor  pressures  above  the  liquid. 

On  the  average,  the  calculated  vapor  pressures  between  238°  and  298°  K  were 
30  percent  lower  than  the  experimental  ones  (fig.  3).  This  difference  is  not  significant, 

however,  inasmuch  as  at  the  low  pressure  levels  involved,  the  average  difference  is 

-  9 

only  in  the  order  of  10  atmosphere.  Furthermore,  the  lower  calculated  vapor  pres- 
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i 


sures  are  consistent  with  the  assumption  of  reference  26  that  the  experimental  vapor 
pressures  are  too  high. 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  April  26,  1965, 
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APPENDIX  -  SYMBOLS 


Ft  -  riQ 


AH° 


rotational  constant,  cm 

heat  capacity  at  constant 
pressure  for  standard 
state,  cal/(mole)(°K) 

equilibrium  specific  heat  at 
constant  pressure  of  mix¬ 
ture,  cal/g 

frozen  specific  heat  at  con¬ 
stant  pressure  of  mixture, 
cal/g 

reactive  specific  heat  at  con¬ 
stant  pressure  of  mixture, 
cal/g 

spectroscopic  constant  for 
rotational  stretching,  cm”^ 

dissociation  energy  at  0*^  K, 
cal/mole 

Gibbs  free  energy  for  stand¬ 
ard  state,  cal/mole 

sensible  free  energy  for 
standard  state,  cal/mole 

sum  of  sensible  enthalpy  at 
T°  K  and  chemical  energy 
at  0°  K  for  standard  state, 
cal/  mole 

sensible  enthalpy  for  stand¬ 
ard  state,  cal/mole 

enthalpy  change  for  forma¬ 
tion  of  substance  from  ele¬ 
ment  in  atomic  gas  state, 
cal/  mole 


(AH^) 

^  f 


h 

^Nrap 


K 


Kf 


M 


n 

P 

Pi 

R 


enthalpy  change  for  formation 
of  substance  from  assigned 
reference  element,  cal/mole 

chemical  energy  at  0°  K  for 
standard  state,  cal/mole 

enthalpy,  cal/g 

heat  of  vaporization  to  1  gram 
of  equilibrium  vapor 

equilibrium  constant  for  re¬ 
action  of  formation  from 
element  in  atomic  gas  state 

equilibrium  constant  for  re¬ 
action  of  formation  from 
assigned  reference  element 

molecular  weight  based  on 
chemical  scale  of  natural 
oxygen,  g/mole 

mole  fraction  of  r”  species, 
moles  of  i  species  per 
mole  of  mixture 

number  of  moles  per  gram  of 
equilibrium  mixture 

total  vapor  pressure,  atm 

th 

partial  pressure  of  r“ 
species,  atm 

universal  gas  constant, 

1.  98726  cal/(mole)(°K)  or 
82.  05971  (cc)(atm)/ 
(mole)(°K) 

entropy  for  standard  state, 
cal/(mole)(°K) 
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s  entropy,  cal/(g)(°K) 

ASyap  entropy  of  vaporization  to 

1  gram  of  equilibrium  vapor 

T  absolute  temperature,  °K 

V  specific  volume  of  vapor,  cc/g 

f  Vi 

x-  weight  fraction  of  species, 

1  XT* 

g  of  species  per  g  of 
mixture 

vibration-rotation  interaction 
e 

constant  for  diatomic  mole¬ 
cule,  cm'^ 

p  density  of  vapor,  g/lOOO  cc 

a  standard  deviation 


CO 


e 


Cl)  x„,co  v„ 
e  e’  e-'e 

Subscripts: 

c 

f 

g 

1 

2 


zero-order  vibrational  fre¬ 
quency  for  diatomic 
molecule,  cm”^ 

anharmonicity  constants  for 
diatomic  molecule,  cm"^ 

condensed  phase  property 

formation  from  assigned 
reference  element 

property  of  equilibrium 
vapor  mixture 

monomer  property 

dimer  property 
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conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
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shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
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